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ABSTRACT 

Hypervelocity stars (HVSs) discovered in the Milky Way (MW) halo are thought to be 
ejected from near the massive black hole (MBH) at the galactic centre. In this paper 
we investigate the spatial and velocity distributions of the HVSs which are expected 
to be similarly produced in the Andromeda galaxy (M31). We consider three different 
HVS production mechanisms: (i) the disruption of stellar binaries by the galactocentric 
MBH; (ii) the ejection of stars by an in-spiraling intermediate mass black hole; and 
(iii) the scattering of stars off a cluster of stellar-mass black holes orbiting around the 
MBH. While the first two mechanisms would produce large numbers of HVSs in M31, 
we show that the third mechanism would not be effective in M31. We numerically 
calculate 1.2 x 10^ trajectories of HVSs from M31 within a simple model of the Local 
Group and hence infer the current distribution of these stars. Gravitational focusing 
of the HVSs by the MW and the diffuse Local Group medium leads to high densities of 
low mass IM©) M31 HVSs near the MW. Within the virialized MW halo, we expect 
there to be of order 1000 HVSs for the first mechanism and a few hundred HVSs for 
the second mechanism; many of these stars should have distinctively large approach 
velocities (< —500 km s~^). In addition, we predict ~ 5 hypervelocity RGB stars 
within the M31 halo which could be identified observationally. Future MW astrometric 
surveys or searches for distant giants could thus find HVSs from M31. 

Key words: stellar dynamics - black hole physics - galaxies: individual: M31 - Local 
Group 



1 INTRODUCTION 

Hypervelocity stars (HVSs) travel at such high speeds that 
they are not gravitationally bound to the galaxy from which 
they originate. The first HVS was discovered in the Milky 
Way (MW) halo by Brown et al. (2005); nine other HVSs 
in the MW halo have since been found by Edelmann et 
al. (2005), Hirsch et al. (2005) and Brown et al. (2006a,b, 
2007b). By selection, most of the observed stars are B-type 
stars. The HVSs have typical speeds of 550 km s~^. There 
is also evidence for a bound population of stars similar to 
the HVSs but with speeds slightly lower than the escape ve- 
locity of the MW (Brown et al. 2007a). The observations of 
MW HVSs indicate that all but one of these stars originate 
from the centre of the galaxjQ. This is easily understood if 
HVSs are ejected by the massive black holes (MBHs) which 
inhabit the centres of galaxy bulges. 

The ejection of HVSs from galactocentric MBHs was 
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^ The HVS discovered by Edelmann et al. (2005) appears to orig- 
inate from the Large Magellanic Cloud. 



first predicted by Hills (1988). The proposed mechanism 
was the disruption of a tight stellar binary by the MBH, 
which results in the capture of one of the stars and the 
ejection of the other as a HVS. This tidal breakup ("TB") 
mechanism has been further investigated by Yu & Tremaine 
(2003), Gualandris et al. (2005), Ginsburg & Loeb (2006), 
and Bromley et al. (2006). Another model for the produc- 
tion of HVSs is the scattering of stars bound to the galactic 
MBH by an intermediate or high mass in-spiraling black 
hole ("IBH" mechanism - Yu & Tremaine 2003; Levin 2006; 
Baumgardt et al. 2006; Sesana et al. 2006, 2007b). Though 
the results of Sesana et al. (2007a) indicate that this second 
model may be inconsistent with the statistics of the observed 
MW HVSs, it cannot be ruled out in other galaxies. The 
third mechanism for the ejection of HVSs was proposed by 
O'Leary & Loeb (2006); this involves the scattering of stars 
by a cluster of orbiting stellar-mass black holes which have 
segregated around the galactocentric MBH ("BHC" mech- 
anism). It has not yet been conclusively determined which 
of these mechanisms is dominant in the MW, and so we will 
consider them all in our work. 

The HVSs originating from the MW centre have been 
studied in some detail by both observers and theorists. How- 



2 B. D. Sherwin, A. Loeb and R. M. O^Leary 



ever, the production of HVSs is of course not unique to our 
galaxy; indeed it is expected to be a feature of all galax- 
ies with central MBHs. The HVSs which escape from other 
galaxies will become, to quote Hills (1988), "intergalactic 
tramps" . In particular, the neighbouring Andromeda galaxy 
(M31) is fairly similar to the MW and has a MBH at its cen- 
tre (Bender et al. 2005). We thus expect a HVS population 
similar to that observed in the MW to be produced in M31, 
with some of the escaping stars gravitationally focused into 
the MW. Observations of such stars could be used to deter- 
mine the HVS production mechanism, elucidate the prop- 
erties and merger history of the core of M31, and constrain 
the distribution of mass within the Local Group. This leads 
to two questions which are the subject of our paper: (i) what 
is the distribution of stars ejected from M31?; and (ii) can 
we observe them from our galaxy? 

To answer these questions one needs to consider the his- 
tory of the MW and M31 within the Local Group of galax- 
ies (of which they are the largest constituents). Kahn and 
Woltjer (1959) first introduced the "timing argument", by 
which the MW and M31 are assumed to have formed close 
to each other in the early universe and were subsequently 
pulled apart by the general cosmological expansion. The two 
galaxies are now decoupled from the expansion and have 
traced out much of a Keplerian orbit, which explains the 
current observed approach velocity of ^ 120 km s~^ (Bin- 
ney & Tremaine 1987). Given the current approach veloc- 
ity, the separation of the galaxies, the age of the universe, 
and assuming that the Local Group has no angular momen- 
tum, one can solve Kepler's equations to find an estimate 
of the Local Group mass - ~ (3-5) X 10^2 M© (see Binney 
& Tremaine 1987, Fich & Tremaine 1991). The timing ar- 
gument has since also been applied to systems with angular 
momentum and realistic mass distributions (see e.g. Peebles 
et al. 1989; Valtonen et al. 1993; Peebles 1994; Peebles et al. 
2001; Sawa & Fujimoto 2005; Loeb et al. 2005). 

The timing argument also allows an approximation to 
be made about the separation of the MW and M31 for the 
dynamics of HVSs. A feature of the galaxies' eccentric or- 
bit is that the galaxies spend a large fraction of the orbital 
time near apogalacticon, with a relative speed less than the 
current value. Since the HVSs move at speeds which are at 
least a few times larger than the current relative radial speed 
of M31 and the MW, a reasonable approximation in follow- 
ing the HVS dynamics is to assume a static configuration 
over the last 10^° years, with M31 and the MW separated 
by their current distance. However, we also verified that the 
HVS distributions are not significantly affected by the ad- 
dition of a radial component to the HVS velocities of order 
100 km s"^ 

Recent studies of the Local Group investigated the 
galactocentric transverse velocity of M31. Peebles et al. 
(2001) applied the action principle to the motion of galaxies 
within 20 Mpc of the Local Group and obtained a trans- 
verse velocity of ~ 200 km s~^. Loeb et al. (2005) analysed 
the proper motion of M31's satellite galaxy M33 using nu- 
merical simulations and estimated the transverse velocity of 
M31 as ^ 100 km s~^. A recent analysis of the motion of 
M31's satellite galaxies by van der Marel & Guhathakurta 
(2007) gave an M31 transverse velocity of ^ 40 km s~^.In 
our work we neglect the transverse velocity as it has only a 
small infiuence on the trajectories of the fast HVSs. 



In the following sections we examine the trajectories of 
HVSs from M31 within the Local Group and estimate their 
current distribution and observability. We first construct a 
simple model of the Local Group and its gravitational field 
as described in §2.1. We then numerically integrate the equa- 
tions of motion for a large number of possible HVS trajecto- 
ries {'^ 10^) to find the current HVS distribution within the 
Local Group, sampling from assumed probability distribu- 
tions of HVS ejection direction, ejection speed, and stellar 
mass, as detailed in §2.2. An analysis of how to perform 
the simulations most efficiently and a discussion of the nu- 
merical methods we use to find the HVS trajectories and 
distributions is given in §3. The results of our simulations 
are presented in §4 and discussed in §5. We comment on 
the observational outlook for our results in §6 and finally 
conclude in §7. 



2 MODELLING THE SYSTEM 
2.1 The Local Group 

The trajectories and the spatial distributions of the HVSs 
produced in M31 depend on the gravitational field within the 
Local Group. A number of models have been constructed of 
the mass distribution within the entire Local Group (Cox 
& Loeb 2007) and of the MW and M31 (Klypin et al. 
2002; Widrow k Dubinski 2005; Seigar et al. 2006). Much 
of the Local Group's mass is contained within the MW and 
M31 halos. These galaxies are surrounded by a diffuse Local 
Group medium of dark matter and gas. Based on cosmo- 
logical simulations in which similar systems form, the mass 
in the diffuse medium is expected be a substantial fraction 
of the total mass (Cox k, Loeb 2007); indeed the timing 
argument's mass estimate of ^ (3-5) x 10^^ is generi- 
cally larger than the sum of the MW & M31 halo masses 
(- 2.6 X lO^^M©). 

We will thus assume in our model that the mass in the 
intragroup medium is equal within a factor / to the mass 
contained within the virialized halos of M31 and the MW, 
Ma + Mm • For simplicity, this diffuse Local Group mass is 
taken to be uniformly distributed within a sphere of radius 
Rg about the midpoint between M31 and the MW. The two 
galaxies are separated by 780 kpc (McConnachie et al. 2005, 
Ribas et al. 2005). 

The mass of both the bulge and disc of M31 and the 
Milky Way is much smaller than the masses of their dark 
matter halos, so the gravitational infiuence of the bulge and 
disc is only important very close to the galaxies' centres. 
As we are interested in approximate HVS trajectories and 
distributions within the entire Local Group and do not re- 
quire precise trajectories near the galactic centre, we neglect 
the gravitational field of the disc and bulge in our Local 
Group model. However, the M31 disc and bulge significantly 
slow an HVS which is escaping the galactic centre; we take 
this into account by reducing the HVS initial velocity by an 
amount which corresponds to the kinetic energy required to 
escape the disc and bulge. To estimate this reduction in ve- 
locity, we assume a Hernquist bulge (Hernquist 1990) with 
a mass ruh = 2 x 10^° M© (Klypin et al. 2002) and a scale 
radius ah = 1.8 kpc (Widrow & Dubinski 2005), as well as 
an exponential disc with a mass rud = S x 10^° and a 
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scale radius Rd — 5.5 kpc (Widrow & Dubinski 2005). The 
central potential of the bulge is —Gmb/at (Widrow & Du- 
binski 2005), and we calculate the central exponential disc 
potential to be —Grrid/ Rd- The escape velocity Vhd from the 
bulge and disk is thus Vhd ~ 470 km s~^; a HVS with an 
initial speed v (once outside the influence of the MBH) has 
a speed vq — \A^^"^-^^ after escaping the bulge and disk. 
We assume that the gravitational field of the M31 disc does 
not modify the initial angular distribution of HVSs as the 
disc's field is small and mainly radial. 

We model the dark matter distribution within the MW 
& M31 galaxy halos with the Navarro, Frenk and White 
(1996; hereafter NFW) density profile which follows from 
cosmological galaxy formation simulations: 



p{r) 



Po 



(1) 



(cr/i?2oo) (l + cr/i?2oo)^' 

where c is the concentration parameter, R200 is the virial 
radius of the galaxy at which we cut off the halos (interior 
to which the average density is 200 times the mean density 
of the Universe), and where po is a constant chosen to give 
the correct total virial mass Ma (for M31) or Mm (for MW). 
The total density in our model is simply the superposition of 
the densities of the two galaxies and the intragroup medium. 

Our Local Group model is thus completely specified by 
the parameters Ma, Mm, /, i?A,2oo, Rm, 200 , Rg, ca and 
Cm- The standard values of these constants which we use 
for our model are Ma = 1.6 x 10^^ M©, Mm 1 x 10^^ M©, 
/ = 1, CA = CM = 12, i?A, 200 - 273 kpc, Rm,2oo - 206 kpc 
and Rg =0.9 Mpc. These numbers are similar to those in 
Cox & Loeb (2007), which are based on observations and are 
also consistent with simulations of the development of the 
Local Group. The model is depicted in Figure [1] In following 
the dynamics of fast HVSs near the turnaround epoch of the 
Local Group, we will assume that the configuration shown 
in Figure [1] has not changed significantly over the age of the 
Local Group (- 10^ Vr). 

The simplicity of our model enables us to find analytic 
expressions for the gravitational field, which increases the 
efficiency of our numerical calculations. Making use of the 
spherical symmetry of all the features in our model, the 
gravitational field can easily be calculated from Gauss' law. 
Within the Local Group the diffuse mass leads to the field 

, ^ -G{MA + MM)f 
gG(r) = -3 r. 

The field within the MW resulting from the NFW halo pro- 
file is a central field of magnitude (Navarro et al. 1996) 

GMm 1 ln(l + cx) - cx/{l + cx) 



|gM(r) 



02 2 

^M,200 ^ 



ln(l + c)-c/(l+c) 



(2) 



where tm is the position of the MW centre and x = |r — 
rM\/RM,2oo- The expression for the M31 field is analogous. 
The fields have a 1/r^ form outside the galaxy or Local 
Group. 

While the halo field is not singular at the centre of the 
galaxy, we avoid numerical difficulties caused by the denom- 
inator by including a softening length cRm, 200'- 

GMm 1 ln(l + cx) -cx/(l + cx) 



r(r) = 



7?2 

^M,200 



(x + e)2 ln(l + c) -c/(l + c) 
- Tm 




|r - fmI + eRA 



(3) 



Figure 1. Our model of the Local Group. The dashed line indi- 
cates the range of our later plots. 



We choose a small value of e = 0.01 which introduces a 
realistic core in the region of the galactic bulge and is found 
not to influence trajectories away from the galaxy centres. 
The total gravitational field is found by superposition: 

gt = gG + gM + gA. 



2.2 The HVS Production Rate 

To accurately model the ejection of HVSs from M31, we 
must determine the number of HVSs produced, the prob- 
ability distribution of the HVSs' initial velocities, and the 
distribution of their stellar masses. For simplicity, we will 
assume that HVSs have been ejected in the same manner 
throughout the lifetime of the Local Group; the velocity and 
mass probability distributions are thus assumed constant 
over the last 10^° years, with no HVS production before this 
time. We neglect the possibility that multiple mechanisms 
simultaneously contribute to the HVS production and only 
consider the mechanisms individually. 

We first discuss the number of HVSs produced in the 
MW. A HVS production rate in the MW can be estimated 
by extrapolating the observations of hypervelocity blue stars 
in the MW halo to lower masses. We assume that star for- 
mation in the MW galactic centre has followed a flat initial 
mass function ^ as found in current observations of the 
galactic centre, |^ oc m~^'^^ for m > IM© (Maness et al. 
2007). The cumulative distribution of stellar masses in the 
galactic centre P(m) is given by P{m) oc (see Alexan- 

der & Sternberg 1999, O'Leary & Loeb 2007, Sesana et al. 
2007a, Perets et al. 2007). This implies that P{m) oc m"^•^^ 
where ti is the stellar lifetime, which is less than a Hubble 
time, and we have used ti oc m~^'^ (Kippenhahn & Weigert 
1990). Taking the masses of ejected HVSs to be distributed 
according to P{m) and assuming that HVSs were produced 
at a constant rate, the HVS sample of Brown et al. (2006a, 
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2007b) leads to a rough estimate of slightly more than 10^ 
HVSs ejected from the MW centre over the age of the local 
group. 

The number of HVSs produced in M31 will not nec- 
essarily be the same as the number produced in the MW. 
While M31 is overall very similar to the MW, the masses 
of the MBHs in the galaxies are very different; the mass 
of the MW's MBH (Sgr A*) is ^ 3.6 x 10^ M© (Eisen- 
hauer et al. 2005), whereas the mass of the MBH in M31 
is ^ 1.4 X 10^ M© (Bender et al. 2005). In addition, the 
stellar populations in the centres of the two galaxies have 
different velocity dispersions, mass functions and densities 
(Bender et al. 2005, Genzel et al. 2003). As the MBH and 
its stellar environment play a large role in all HVS produc- 
tion mechanisms, we must investigate how the differences 
will influence the relative HVS production rates in the MW 
and in M31. We can then simply scale the rates in the MW 
(inferred from observations) for each of the mechanisms to 
obtain the total number of HVSs produced in M31. 

We briefly review some features of the stellar environ- 
ment of MBHs. The stars in the galactic centre have a typ- 
ical mass TTis ~ IMq. The radius of influence, r^, is the 
distance to which the MBH significantly influences the ve- 
locity and spatial distribution of stars; at this radius the 
kinetic energy in the random motions of the stars (char- 
acterised by the velocity dispersion a) is roughly equal to 
the stars' potential energy in the gravitational field of the 
MBH, so that n ~ GM./crf , where M, is the mass of the 
central MBH. Outside the radius of influence, a is approx- 
imately constant with radius {a ^ cr^); inside the radius of 
influence, the velo city disp ersion roughly equals the circular 
velocity (cr(r) ^ GM^jr). The density of stellar mass near 
the black hole is assumed to be described inside the radius 
of influence by the theoretically derived expression (Young 
1980): 



p(r) = pi{r/ri) 



-1.5 



(4) 



Outside the radius of influence, we assume an isothermal 
stellar distribution, so that 

-2 



p(r) = piir/nY 



(5) 



While one could also determine a density profile based solely 
on current observations and not on theoretical considera- 
tions, the current density would not necessarily correctly 
describe the density over the past 10^° years; we believe 
that the above theoretical expressions for the density should 
provide a better average model over the age of the Local 
Group. Our density profile agrees fairly well with observa- 
tions of the galactic centres of both the MW (Genzel et al. 
2003) and M31 (Lauer et al. 1993). We use the following 
values of r^, pi and Gi (assuming a density profile as in Eq. 
SI): for the MW tm = 1.66pc, = 6 x lO^Mo/pc^ and 
CTM ^ 100 km s"^ (Genzel et al. 2003); for M31 ta = 30pc, 
= 6 X lO^Mo/pc^ and a a ~ 150 km s~^ (Bender et al. 
2005, Lauer et al. 1993). We assume a binary fraction 77 of 
ry = 0.1. 

In the BHC model, a cluster of stellar-mass black holes 
has segregated close to the MBH. The stars near the MBH 
are moving at high speeds, so that scattering of these stars 
off the orbiting stellar-mass black holes leads to hyperve- 
locity ejections. The BHC mechanism of course presupposes 
that the stellar-mass black holes have had sufficient time 



since the formation of the galaxy to segregate into a central 
cusp. The segregation occurs by dynamical friction on the 
timescale td ~ —tr, where tr is the relaxation timescale 



tr{r) = 0.34 



G^m,ap{r) In A 



(6) 



(Binney & Tremaine 1987), rrih ^ 7Mq is the mass of a 
typical segregating black hole and In A is the Coulomb log- 
arithm which we take as ~ 10. Using the expressions for 
density and velocity dispersion above, it is easily shown that 
the relaxation and dynamical friction timescales have a con- 
stant minimum value below the radius of influence and rise 
as outside it. Thus we can calculate the minimum time 
required to form a cluster of black holes from the minimum 
dynamical friction timescale at the radius of influence td{ri). 
For the MW, the minimum dynamical friction timescale is 
tM ~ 1.4Gyr (Miralda-Escude k Gould 2000) which implies 
that a black hole cusp would form within a small fraction 
of the age of our galaxy and produce HVSs. Equation JB} 
implies that the minimum dynamical friction timescale in 
M31 is 



tA ^ 



\cns/[ y 



Pm 

PA 



(7) 



which gives Ia ~ 340tM ~ 480Gyr. This greatly exceeds the 
age of the universe, and it follows that no central black hole 
cluster has formed in M3lEl. We thus do not expect any HVS 
production in M31 by the BHC mechanism, and we will not 
consider this mechanism in our calculations. 

We now consider the effect of the increased MBH mass 
and different stellar environment in M31 on the TB mech- 
anism. In the TB mechanism, a binary which approaches 
the MBH to within the tidal radius rt is disrupted; one of 
the stars is captured and its binary companion is ejected 
at high velocity. A binary will be disrupted within an or- 
bital period tp if its orbital specific angular momentum (per 
unit mass) j is sufficiently small, j < jmin ~ \/GM.rt ~ 
^^GM.a{M./msY/^, where a is the binary separation. This 
range corresponds to a region in the angular momentum - 
energy phase space which is known as the "loss-cone". The 
rate of binary disruption is either limited by the orbital pe- 
riod ("full loss-cone") or by the time for binaries to diffuse 
into the loss cone ("empty loss cone"); to calculate the to- 
tal rate one must consider both these regimes (see Frank 
& Rees 1976, Lightman & Shapiro 1977, Cohn & Kulsrud 
1978, Magorrian k Tremaine 1999, Yu k Tremaine 2003 
and Perets et al. 2007). The cross-over from the empty to 
full loss-cone regimes occurs at a critical radius Vc where 
the root-mean-square angular momentum change Aj due to 
stellar encounters during one orbital period tp is equal to 
jmin- We can calculate the critical radius as follows. Vc is 
defined by 



(8) 



^ A caveat to this conclusion involves violent relaxation (Lynden- 
Bell 1967) which could result from a merger event. However, the 
effect of a merger is expected to be episodic and would not lead 
to a cumulative segregation process that is considered here over 
long periods of time. 
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The root mean square specific angular momentum transfer 
over one orbital period is (Lightman & Shapiro 1977) 



/ \ 1/2 



Jc 



(9) 



where jcirc{r) ~ a{r)r is the specific angular momentum of 
a star on a circular orbit. This implies 



_ tp(rc) 

Jcirc (Vc) J tr{rc) ' 



(10) 



Noting that within the radius of influence a is of order the 
circular velocity and that outside the radius of influence it 
is a constant, we take tp ^ 27vr/a{r). Using the expressions 
given previously for jmin, jcirc and tr, we obtain an equation 
involving rc 

^1/3 

(11) 



GM.a{M./ms 



2TircG'^msp(rc 



0.034cr4 



Solving this equation using a = 0.1 AU leads to a critical 
radius significantly larger than the radius of influence for 
both M31 and the MW {n/rc < 1). 

We now make an order of magnitude estimate of the 
HVS production rates using our density profile. As is dis- 
cussed in Perets et al. (2007) §4, the contribution to the 
total rate V is dF/dr oc r^^^ for r < ri, dV/dr oc for 
Ti < r < Tc and dF/dr oc for r > Vc- The continuity 
of dF/dr allows us to estimate the relative contributions of 
these three different regions to the total rate; it is easily 
shown that the contribution from the full loss-cone region 
r > rc is smaller than that from the empty loss cone region 
r < rc by a factor ri/vc and is thus negligible. Perets et 
al. (2007) give an approximate expression for dF/dr in the 
empty loss cone regime (their Eq. 17): 



(12) 



dFc ^ Ns{<r) 
dr rtr(rc) 

where Ns(< r) is the number of binaries enclosed within r. 
The total rate can thus be obtained from an integration over 
r: 

407rG2r?p?rf 



0.034 X 6(7? 



(13) 



Equation ([13]) allows us to find the expected ratio of HVS 
production in the MW and M31 based on the observational 
data mentioned abovjf]: 



Fm \PM J xtm) \cfa. 



' 0.2 



(14) 



This rough estimate indicates that the TB mechanism pro- 
duces 1-10 times fewer HVSs in M31 than in the MW. 

The IBH mechanism is also affected by the differences 
between the MW and M31. The HVS ejection rate was cal- 
culated for M31 by Lu et al. (2007) by adopting the method 
of Yu & Tremaine (2003); assuming an IBH/MBH mass ra- 
tio of = 0.01 and an IBH-MBH axis of 5.9mpc the rate is 
approximately an order of magnitude lower than in the MW 
(of order 10" Vy^)- 

^ We also verified that this simple formula gives a rate very sim- 
ilar to the more exact Yu & Tremaine 2003 calculation for the 
MW. 



To summarise, the BHC mechanism produces no HVSs 
in M31, and both the TB and IBH HVS production rates 
are suppressed in M31 by a factor of order 1-10 compared 
to the MW. However, this may be partially compensated for 
by the somewhat higher ejected star velocities in M31 (Yu 
& Tremaine 2003) which decrease the proportion of bound 
ejected stars. The number of HVSs produced over the past 
10^° years in the MW is expected to be somewhat above 
10^. For our calculations with both these mechanisms we 
thus assume that the total number of HVSs produced in 
M31 is 10^. This estimate is fairly uncertain and relies on 
the assumptions that the stars are isotropically distributed 
and that the effects of any massive perturbers (Perets et 
al. 2007) are equal in the MW and in M31. The number of 
HVSs produced could easily be an order of magnitude lower 
or higher than our estimate. However, our statistical results 
can be simply scaled to a different number of ejected stars 
based on future observations and calculations. 



2.3 HVS Velocities and Masses 

We now consider the probability distributions for the HVSs' 
ejection velocities and masses. Both the TB and IBH models 
have the welcome feature that the form of the probability 
distribution for the velocity of an ejected star is at most 
weakly dependent on its mass: for the IBH mechanism, the 
velocity is entirely independent of the mass of the HVS (as 
the HVS acts as a test particle in the gravitational field of 
the much heavier black holes); for the TB mechanism, the 
velocity of a HVS of mass mi shows a very weak dependence 
oc y/rn2{mi + 7712)""^^^, where m2 is the mass of the other 
binary component (Hills 1988). We thus approximate the 
velocity and mass distributions as independent, which makes 
our numerical exploration of their statistics much simpler. 

For the TB model, Sesana et al. (2007a) apply the re- 
sults of Bromley et al. (2006) to binaries with a semi-major 
axis a distributed uniformly in In a (Heacox 1998). Following 
their Figure 1 (TBf model), the TB mechanism gives rise to 
a distribution of star ejection velocities v (after escaping the 
gravitational field of the MBH) of approximately 



P{v) 



0, 



V > Vl 

V < Vl 



(15) 



where P{v)dv provides the probability that a star has an 
ejection speed in the interval {v,v + dv). While this veloc- 
ity distribution was derived for the MW, we will assume 
that the power law index of the distribution is characteris- 
tic of the HVS production mechanism and does not depend 
strongly on the MBH mass, so that this distribution can also 
be applied to M31. vi is the lower cutoff velocity, which we 
define as the maximum speed with which stars do not escape 
the M31 halo over 10^° years for any ejection direction, vi 
is thus very similar to (but slightly smaller than) the escape 
velocity. 

The HVS velocity distribution for the IBH model was 
calculated numerically for the MW by Sesana et al. (2006, 
2007a,b). With an IBH/MBH mass ratio of z/ = 1/729 and 
an initial orbital eccentricity of e = 0.9, they consider the 
entire in-spiral of the black hole and obtain a HVS velocity 
distribution oc v~^'^. This distribution, with the high ve- 
locity tail, is nearly independent of the assumed values of 
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v and e (A. Sesana, private communication). Again, we as- 
sume that the M31 HVS ejection velocity distribution has 
the same power law index as was found in the MW and thus 
obtain for M31 



P{v) oc 



0, 



V > Vl 

V < Vl 



(16) 



As explained previously, the initial speed vq we us e for our 
simulations of HVS trajectories is given by vq — — v^^ 
(taking into account the slowing by the M31 bulge and disc.) 

The dependence of HVS ejection probability on stellar 
mass is very weak for both the TB and IBH mechanisms, 
and so we will assume that the distribution of ejected stel- 
lar masses follows the distribution of main sequence stars 
in the galactic centre. However, the star formation history 
in M31 over the past 10^*^ years is not known, and so for 
simplicity we will assume a Salpeter initial mass function 
4^ oc m~^'^^ . As observations indicate that the formation 

dm 

of low mass stars is suppressed in galactic centres (Maness 
et al. 2007) and very low mass stars are difficult to observe, 
we cut off the mass distribution below IMq. However, since 
the lifetime of IMq stars is roughly the same as the Local 
Group age, the distribution of the longer lived lower mass 
stars can be obtained by simply rescaling the IMq distribu- 
tion to reflect any differences in the total number of stars. 
As discussed previously, taking stellar lifetime into consid- 
eration implies that the HVSs ejected from M31 thus have 
masses distributed according to 



P{m) oc 



m-^•^^ m>lM0 
0, m < IMq 



(17) 



We assume that the HVSs are ejected isotropically from the 
M31 centre. 



3 NUMERICAL METHODS AND ANALYSIS 
OF OUR MODEL 

3.1 Numerical Methods 

In order to obtain HVS distributions and velocities we must 
simulate a large number of individual stellar trajectories. 
The basic challenge of this calculation is to integrate New- 
ton's gravitational force law. The integration method we use 
is the Verlet method (Verlet 1967), which has a number of 
advantages: it is symmetric (i.e. it uses information from 
both ends of a step), it is a fourth order method, and it is 
sufficiently fast to allow a large number of trajectories to 
be calculated; many HVS trajectories are needed to obtain 
statistically valid information. A fixed step size is used for 
our integration as this greatly simplifies the determination 
of the HVS distributions. 

To find one HVS trajectory, our integrator was run for a 
time corresponding to the age of the Local Group, which we 
take to be 10^*^ years. We tested our program by comparing 
its results with a number of analytically calculated trajec- 
tories for different force laws. Our numerical results agreed 
with the analytical solutions for all the configurations we 
tested (e.g. Keplerian orbits, simple harmonic motion). For 
our numerical simulations, we chose to use 10^ integration 
steps per trajectory. This number of steps was not compu- 
tationally expensive and gave results which agreed very well 



with analytic solutions. It was verified that the trajectories 
were identical for a wide range of step-numbers (from 500 
to 2 X 10^), which suggests that the integration scheme is 
stable and that our results are reliable. 



3.2 Symmetry Considerations 

Our model for the Local Group has a number of important 
symmetries, which can be exploited to make the necessary 
numerical calculations more efficient. As is seen in Figure 
[ll the system exhibits rotational symmetry along the z axis 
connecting the two galaxies, which means that every tra- 
jectory of an HVS from the centre of M31 follows a path 
which lies in a single plane through the inter-galactic axis. 
This feature allows us to reduce our problem from three to 
two dimensions (and thereby increase the computational ef- 
ficiency), by "folding" all possible trajectories into one cross- 
section plane within which we perform all our calculations. 
To regain three dimensional results we can simply deproject 
our values. An isotropic ejection in three dimensions corre- 
sponds to the following distribution in two dimensions: 



P{e) = -smt 



(18) 



where 6 is the angle between the intergalactic 2;- axis and the 
initial HVS velocity vector. 

In addition, the reflection symmetry of the upper and 
lower half planes means that a "reflected" trajectory is also 
a valid one. We can thus limit ourselves to sampling trajec- 
tories that begin in the upper half plane. 

The final and most important symmetry involves time. 
In our model the mass distribution as shown in Figure [1] is 
static. This means that every point (not just the endpoint) 
of a trajectory is a possible current location of a HVS. For 
every simulated trajectory we thus obtain not one, but many 
possible current star positions. 

The choice of a two-dimensional static Local Group 
model over a three-dimensional dynamic model reduces the 
computational time for the HVS distributions by a factor of 

10^. 



3.3 Calculation of the HVS Distributions 

To find the HVS distributions and velocities, we must calcu- 
late a large number of HVS trajectories. The initial condi- 
tions for each trajectory are determined by sampling ran- 
domly from the distributions for the HVS ejection angle 
(Eq. [HI) and speed (Eq. [T5p6l). We use Park and Miller's 
(1988) "Minimal Standard" random number generator to- 
gether with the transformation method detailed in Press et 
al. (1996) to obtain initial conditions with the required dis- 
tribution. 1.2 X 10^ trajectories were thus calculated, which 
correspond to 1100 independent samples in both angle and 
velocity. For each trajectory we used 10^ timesteps. This 
means that we obtain a total of 1.2 x 10^*^ points, each of 
which is a possible current position of a HVS. However, due 
to the finite stellar main-sequence lifetime, each of these 
points does not carry equal weight when determining the 
distribution of HVSs within the Local Group. 

We now consider how to obtain the distribution of HVSs 
from the raw trajectory data. The lifetime ti of a star of mass 
m is assumed to be ti — lOGyr(m/M0)~^"^ (Kippenhahn & 
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Weigert 1990). We take the age of the Local Group to be 
^age = 10^° yr. Every point on a trajectory corresponds to an 
H VS travel time tt which can range from zero to tage • Before 
being ejected from the centre of M31, a star is assumed to 
have resided in the galactic centre for a "dwell time" td- For a 
star with mass m and corresponding main-sequence lifetime 
tl, a point along a trajectory corresponding to a travel time 
tt can contain a "live" star if tt -\- td < ti. Therefore the 
probability F{tt) of a star at a point corresponding to tt 
being "live" is equal to the probability that td < {ti — tt). 

To determine this probability we must consider the dis- 
tribution of dwell times P{td). There are two conditions on 
td. A "live" star cannot have resided in the galactic centre 
for more than its main-sequence lifetime t^, so that td <ti. 
Of course, it can also not have been there before the for- 
mation of the Local Group tage = 10^° years ago, which 
implies td + tt < tage- Thus the maximum value of td is 
tcZ,max =min(ti,tage — tt) . Wc assumc that a live star is 
equally likely to have been in the galactic centre for any 
time up to the maximum permissible value tcz,max, and the 
dwell time thus obeys a uniform probability distribution 
P(td) = l/td,max between and td,max. 

The probability F(tt) of a star at a point with a travel 
time tt being live is equal to the likelihood of fulfilling the 
condition td < {ti — tt), i.e. 

F{tt)= r P{td)dtd (19) 

^0 

The different ranges of this integral imply different 
forms for P{td), which complicates this (otherwise straight- 
forward) integral. If ti < 0. Stage 

If tl > 0.5iage 

r tt<te.ge-tl 
F(tt) = <^ 1 - tl>tt> tage - tl (21) 

t 0, tt> tl 

The density of HVSs within the Local Group is thus 
calculated as follows. We partition the ejected stars into 
six mass bins containing equal numbers of HVSs according 
to the mass distribution given by equation ()17|) . The mass 
ranges corresponding to these bins are (1 — l.O5)M0, (1.05 — 
1.11)M0, (1.11-1.2O)M0, (1.2O-1.33)M0, (1.33- 1.59)M0 
and (> 1.59)M0 for mass bins one through six. The average 
mass in each bin is used to determine the corresponding ti, 
which allows the calculation of F(tt) for all values of tt. We 
then divide the upper half plane of Figure [T] into an array of 
(300 X 600) cells and add up the number of points in each 
cell, with each point weighted by F{tt). Dividing by the to- 
tal number of trajectory points (1.2 x 10^ x 10^) as well as 
the volume corresponding to each cell and multiplying by 
the total number of HVSs produced gives a number density 
n(d,z) of stars/kpc^ for each mass bin; the coordinates (d, 
z) are defined in Figure [1] 

In addition to the distributions for the six mass bins, we 
also consider massive stars with masses above 3M0. These 
are important observationally due to their high luminosity 
(the observed MW HVSs have mainly been 3 — AMq stars). 



We take their lifetime to be represented by ti = 0.35Gyr 
(Schaller et al. 1992). 



4 RESULTS 

4.1 Star Trajectories 

We first investigated the shapes and characteristics of indi- 
vidual star trajectories. Single HVSs were expelled from the 
centre of M31 with a variety of initial conditions and their 
subsequent paths were found by integration of the gravita- 
tional force. For large ejection velocities {v > 900 km s~^) 
the stars escaped the Local Group entirely, without being de- 
flected significantly from their original straight trajectories. 
For velocities smaller than the lower cutoff velocity, the stars 
were bound to M31 and precessed about its centre in loop 
orbits. In addition to the bound and expelled stars, stars 
released under certain angles with a small range of veloci- 
ties ~ 650 — 900 km s~^ were gravitationally pulled into the 
MW or towards the intergalactic axis. Their trajectories had 
very low average speeds of approximately 100 — 200 km s~^, 
as these stars have just enough energy to escape M31's at- 
traction and thus spend a large amount of time at a low 
velocity. 

The value of the lower cutoff velocity was determined as 
follows. We ejected a large number of stars out of the centre 
of M31 with a fixed speed but in different directions, and 
observed whether or not all stars remained within the M31 
halo over the age of the Local Group. Stars released towards 
the centre of the Local Group escaped M31 with lower ve- 
locities than ones released away from the centre. By trial 
and error, we were able to determine that the highest speed 
at which all stars remained bound was vi = 650 km s~^. 

An interesting feature of the system emerged during our 
analysis. It was found that the shape of the trajectories was 
very sensitive to the initial conditions for some parameter 
values. In some cases, a change of less than one percent in the 
starting velocities lead to the ejected stars following entirely 
different trajectories. While the paths with slightly different 
starting conditions were identical initially, it was found that 
the trajectories split after a certain time. These character- 
istics are indicative of chaotic behaviour in some regions of 
phase space. It was found that such chaotic behaviour was 
most common for trajectories which passed through the re- 
gion between the two galaxies. The chaotic sensitivity to 
boundary conditions appeared to require the presence of 
three gravitating objects (the two galaxies and the diffuse 
group medium); when only two of these objects were present, 
the system no longer showed chaotic behaviour. 

While it is difficult to make precise calculations of the 
trajectories of individual stars in a chaotic system, our nu- 
merical experiments demonstrated that statistical proper- 
ties, such as HVS spatial and velocity distributions, were 
not sensitive to slight changes in the initial conditions. 



4.2 HVS Spatial Distributions 

To determine whether the HVSs from M31 are observable, 
it is necessary to find the HVS distribution within the Local 
Group. By integrating 1.2 x 10^ trajectories numerically and 



8 B. D. Sherwin, A. Loeb and R. M. O'Leary 




(a) 



of the number of stars per cubic kpc (TB) 





(b) logiQ of the number of stars per cubic kpc (IBH) 




-200 200 
z/kpc 



1000 
800 
- 600 
400 
200 


-200 
-400 
-600 
-800 
-1000 



(c) average radial velocities in km s ^ (TB) 



(d) average radial velocities in km s ^ (IBH) 




-1000 -600 -200 200 600 1000 



(e) average transverse velocities in km s ^ (TB) 




(f) average transverse velocities in km s ^ (IBH) 



Figure 2. Main-sequence HVS number densities, radial velocities and transverse velocities (with respect to the MW galactic centre) 
for stellar masses of (1-1.O5)M0. Plots are shown for the TB mechanism {Left column) and the IBH mechanism (Right column). White 
colour indicates a value off the scale on the bar to the right of the plots. The positions of the MW and M31 centres are marked with 
black dots, with the MW on the left. See Figure 1 for the definition of the (d, z) coordinates. 
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processing the points on the trajectory as described previ- 
ously, the HVS distributions were found for both the TB 
and IBH mechanisms. We obtain distributions and average 
radial and transverse velocities for all six mass bins and the 
massive stars; Figures 2 and 3 show the results for the first 
and third mass bins respectively. Distributions and velocities 
for the (> 3Mq) stars are depicted in Figure 4. 

The plots of the HVS distributions in Figures 2 and 3 
exhibit a number of interesting features. We first consider 
the TB mechanism mass bins (left columns of Figs. 2, 3). It is 
clear that especially for the lower mass stars, the HVS num- 
ber density departs significantly from an isotropic distribu- 
tion. For the low mass stars of the first bin, there is a sharp- 
featured arc of increased number density (n ^ lO^^/kpc"^, 
roughly 10 times the background density) which extends 
from M31 to the region behind the MW. In addition, the star 
density is very high on axis 10~^/kpc^). For higher mass 
stars, the deviations from isotropy are smaller: for the third 
mass bin, the arc of raised probability is hardly visible and 
the axial maximum is less pronounced; for the fifth and sixth 
mass bins, the HVS distribution is almost entirely isotropic. 
Another feature which is (barely) visible in these plots is the 
presence of shell- like peaks in density around M31. While 
shell-like features are often seen in galaxy merger simula- 
tions (e.g. Hernquist & Quinn 1989), changing our velocity 
distribution showed that these shells were an artificial fea- 
ture due to the sharp cutoff in the lower ejection velocity. 
The density of the > 3Mq stars (Fig. 4) was found to fall 
off quickly and isotropically with distance from M31. The 
density is 10~^/kpc^ at a distance of « 200kpc from M31, 
and only 10~^/kpc^ at a distance of ~ 500kpc. 

The lower-mass HVS distributions for the IBH mecha- 
nism (right columns of Figs. 2, 3) are similar to the corre- 
sponding plots for the TB mechanism, with only small dif- 
ferences. It can be seen that the arc of first-mass-bin stars is 
less pronounced for the IBH mechanism, and that the IBH 
mechanism produces a somewhat lower HVS density in the 
MW halo. The distribution of the > 3Mq IBH stars is sig- 
nificantly different from the TB mechanism distribution. For 
the IBH mechanism, the density of these stars falls off much 
more slowly with distance to give a value of lO^^/kpc*^ at a 
distance from M31 of ~ SOOkpc. 

4.3 HVS Velocities 

We now consider the stellar velocities measured relative to 
the MW centre. It must be noted when considering M31 
HVSs close to the Earth that we have neglected the MW 
bulge and disc, so that the actual M31 HVS velocities there 
will be somewhat higher than we depict in our plots. 

We first describe our results for the TB mechanism. 
Figure 2 shows that for the first mass bin the radial and 
transverse speeds near the MW are typically ~ 400 km s~^. 
It is noticeable that the typical speeds of higher mass stars 
near the MW were substantially greater. For the third mass 
bin (Fig. 3) the typical velocity is ^ 500 km s~^ near the 
MW; for the fifth mass bin it is ~ 700 km s~^ . The high mass 
stars in Figure 4 appear to have particularly high velocities 
of ^ 3000 km s~^ near the MW, but of course the density of 
such stars is negligible in the regions where the speeds are 
so high. 

The velocities resulting from the IBH mechanism (Figs. 



2,3; right column) were significantly higher than for the TB 
mechanism, with typical speeds near the MW of 600 km s~^ , 
800 km s~^, 1000 km s~^ for the first, third and fifth mass 
bins respectively. 

In all velocity plots the average velocities were larger 
at greater distances from the centre of M31 (though the 
numbers of HVSs there were smaller). The highest radial 
speeds occur near the intergalactic axis, whereas the highest 
transverse speeds can be found perpendicular to this axis as 
expected. 



5 DISCUSSION 

In our investigation, we found a chaotic sensitivity of the 
HVS trajectories to initial conditions. Chaotic regions of 
phase space are common for the motion of a test particle 
in a complicated three centre potential; chaos can even be 
found in non- stationary simple two body potentials, such as 
binary systems (see Holman & Wiegert 1999, Musielak et 
al. 2005). 

The arc- like regions of raised density, which are seen 
for low mass stars with both ejection mechanisms, are due 
to the focusing of ejected HVSs by the MW's gravitational 
field. The stars that are focused in this way by the MW have 
only just enough kinetic energy to escape the attraction of 
M31. They are thus very slow and spend a large amount 
of time in the observed "arc", which leads to an increased 
density there, before falling back towards the intergalactic 
axis. While one expects regions of increased density, it is 
unclear why there is a sharp- featured substructure within 
the arc, as there are no sharp steps in the potential near this 
region. Such arc-like features are not observed for higher- 
mass stars. This follows from the reduced main-sequence 
lifetime of these stars (ti oc m~^'^); high mass stars simply 
do not survive long enough to be visible along such long- 
lasting trajectories. The arcs are less pronounced for the 
IBH mechanism because it ejects HVSs at higher speeds, as 
can be seen from Equation ([16]). 

It must be noted that the approximation of a static 
Local Group is not applicable to the longest-lasting trajec- 
tories. Taking the small initial separation of M31 and the 
MW into account would increase the number of low mass 
HVS from M31 we expect to find in the MW halo; the ac- 
tual number densities of M31 HVSs in the MW may thus be 
higher than we calculate in Figure 2. 

A large overdensity of stars was found near the inter- 
galactic z-axis, which is also due to gravitational focusing. 
As described before, the gravity of the MW and the diffuse 
Local Group mass draw HVSs towards the intergalactic axis; 
the axial symmetry of the entire configuration means that 
every such trajectory eventually crosses the axis at some 
point, which increases the density there. It is important to 
note that this gravitational focusing has a very large effect; 
in Figure 2 it can be seen that it increases the number den- 
sity of stars near the the MW by roughly three orders of 
magnitude compared with the number of HVSs at an equal 
M31-distance off axis. 

Of course, the real Local Group mass distribution will 
not be exactly axisymmetric. For instance, the symmetry of 
our model could be broken by the galactic discs, any triaxi- 
ality of the MW or M31 halo (see Gnedin et al. 2005, Yu & 
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(e) average transverse velocities in km s ^ (TB) 



(f) average transverse velocities in km s ^ (IBH) 



Figure 3. Same as Fig. 2 for stellar masses of (1.11-1.2O)M0 . 



Madau 2007), or by other Local Group galaxies. However, as 
the HVSs are usually moving very rapidly, their trajectories 
depend mainly on the large scale mass distribution. While 
local asymmetries may blur the axial overdensity somewhat, 
it is thus unlikely that they would eliminate this feature en- 
tirely. 

The typical velocities of higher mass main-sequence 



stars are larger by selection; as more massive stars have 
shorter lifetimes, they require faster speeds to travel a cer- 
tain distance within their lifetime. Similarly, HVSs have 
larger velocities at greater distances from M31 because only 
a star with a large velocity can reach a location far from 
M31 within its lifetime. The increased isotropy of the distri- 
bution of more massive stars is also due to the finite stellar 
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Figure 4. Same as Fig. 2 

lifetime: more massive live stars are faster, faster stars are 
deflected less by the gravity of the Local Group, and less 
deflection implies a more isotropic distribution. 
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(f) average transverse velocities in km s~-^ (IBH) 
stellar masses above 3Mq. 

from M31 is due to many stars with velocities just below the 
escape velocity turning around at this radius. 



The forms of the radial and transverse velocity distri- 
butions are as expected for stars travelling away from M31. 
The shell-like minimum in the transverse velocity 300 kpc 
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6 OBSERVATIONAL OUTLOOK 

We may now examine the observability of HVSs from M31 
based on the results of our simulations. We first consider the 
large number of low mass (~ IM©) stars expected within 
the MW halo, especially near the intergalactic axis. The 
number density of HVSs in a large part of the MW halo 
is roughly ~ 10"^^"'^'* /kpc^ for the TB mechanism, and a 
few times less for the IBH mechanism. One therefore ex- 
pects of order - 1000 HVSs from M31 in the MW halo for 
the TB mechanism, and a few hundred HVSs for the IBH 
mechanism. These stars have a number of features which 
distinguish them from the background MW stars. First, the 
metallicity of stars from galactic nuclei is different from stars 
in other regions of the galaxy. Second, the HVSs from M31 
have much larger typical speeds than other stars in the MW 
halo. Their average radial speeds near the MW are at least 
400 km s"^ (TB) and 600 km s"^ (IBH); however, some of 
these stars will of course have significantly larger speeds. In 
contrast to HVSs originating in the MW, many M31 HVSs 
will be moving towards the galaxy centre (on the side of 
the MW which faces M31, which is where the solar sys- 
tem is located) with radial speeds significantly exceeding 
the MW escape velocity. However, there is of course a con- 
trast problem in observing such M31 HVSs, as the MW con- 
tains a very large number of ^ IMq stars. A targeted spec- 
troscopic search would probably be unfeasible. However, it 
may be possible to detect these HVSs in planned large scale 
surveys using telescopes such as Pan-STARRS, the Large 
Synoptic Survey Telescope (LSST) or the Global Astromet- 
ric Interferometer for Astrophysics (GAIA) satellite. GAIA 
will measure proper and radial motions of IMq stars out to 
at least lOkpc (Ferryman 2002, Yu k Tremaine 2003) and 
should thus find a small number (of order ^ 10) of M31 
HVSs. LSST will survey roughly half the sky and detect 
IMq stars out to about lOOkpc (Claver et al. 2004, Sesana 
et al. 2007a); it should thus detect a few hundred M31 HVSs 
for the TB mechanism and approximately 100 for the IBH 
mechanism. Even a few stars with sufficiently large nega- 
tive radial velocities or transverse velocities would suggest 
a HVS production in M31. Such observations could be used 
to constrain the characteristics of the M31 galactic centre 
environment as well as the Local Group mass distribution. 
One could perhaps also distinguish between stars from the 
IBH and TB mechanisms by their measured velocities; high 
velocity HVSs produced by the IBH mechanism would in- 
dicate that the M31 MBH has an intermediate mass black 
hole companion (which may be difficult to detect by other 
means) . 

The current population of MW HVSs was discovered 
by observations of luminous main-sequence blue stars in the 
MW halo. From Figure 4, it is clear that the density of 
> 3Mq HVSs from M31 is only significant close to M31 due 
to the stars' short lifetime. Even massive 3 — 8Mq main se- 
quence HVSs would be too faint to be observed near M31 
(unevolved 3Mq blue have an apparent V magnitude of ~ 24 
at the distance of M31). Main sequence HVSs of masses 
greater than ^ 8Mq may be detectable in M31, but one 
would expect only ~ 1 such ejected star to exist at the 
present time (it would be close to the M31 galactic centre). 
An observation of such a HVS would thus be very difficult. 

However, there may be another possibility to detect 
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Figure 5. Number of luminous hypervelocity RGB stars per 
square degree within the distance range of 600-lOOOkpc from the 
MW, as a function of the angle from the centre of M31. Solid line: 
TB mechanism; long-dashed line: IBH mechanism; short-dashed 
line: for comparison, 0.92 x (^/degrees) "-"^ stars per square degree. 
One degree spans a distance of 13.6 kpc near M31. 
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Figure 6. Same as in Fig. [5] but for identifiable hypervelocity 
RGB stars with radial approach velocities < —540 km s"-*^ (see 
text for details on velocity cutoff). Solid line: TB mechanism; 
long-dashed line: IBH mechanism; dotted line: for comparison, 
0.08 X {6/degTees)~^ stars per square degree; short-dashed line: 
0.92 X (0/degrees)-^ as in Fig. [5] 

more distant HVSfl Stars with masses over IMq brighten 
to ~ 10"^ L0 for about 5 x 10^ years when they reach the tip 
of the red-giant-branch (RGB) (Schaller et al. 1992); the lu- 
minosity of these stars would make them observable to very 
large distances. For instance, such a star near M31 has an 
apparent V magnitude of ^ 21.5 and is thus bright enough 
to investigate spectroscopically. We calculated the distribu- 
tion of such giant HVSs (as for the main sequence HVSs) 
and found that they were only present at significant densi- 
ties within ^ 200 kpc of M31 (as such stars are very rare 
due to the short duration of the bright phase, their densities 
are negligible far from their source). 



We are grateful to Nelson Caldwell for this suggestion. 
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To examine the possibility of observing such RGB HVSs 
in the M31 halo, we constructed Figure O which depicts 
the density of hypervelocity red-giant-branch-tip stars per 
square degree as a function of the angular separation from 
M31 (along with an analytic fit to the results). The angu- 
lar density of these stars is seen to fall off as 1/6 at small 
angles; for ^ > 15° the density drops off more rapidly due 
to the finite stellar lifetimes. Significant numbers of RGB 
HVSs can be found at angles beyond the galactic disc and 
bulge (> 1°), with of order ~ 70 (TB mechanism) or ~ 30 
(IBH mechanism) hypervelocity RGB stars in the M31 halo 
at angles of 1° < ^ < 10°. However, to distinguish such 
HVSs from background M31 halo stars and from MW fore- 
ground stars, they need to have radial approach velocities of 
< —420km s~^ (the velocity dispersion of halo stars is less 
than ^ 140 km s~^, see Chapman et al. 2006); it should be 
noted that this cutoff assumes a static Local Group and does 
not include the M31 approach velocity, so that the cutoff ac- 
tually corresponds to a measured galactocentric radial speed 
of roughly < —540km s~^. The angular density of such iden- 
tifiable RGB HVSs is depicted in Figure [6l their density is 
roughly a factor of ten lower, so that one could find of order 
^ 5 such giant stars within the M31 halo (1° < 6* < 10°). 
The background density of non-HVS RGB stars in M31 
(Ibata et al. 2007, Fig. 18) is approximately four orders of 
magnitude greater than that of the HVS giants for ^ > 1°. A 
number of studies of the M31 halo population have recently 
been performed, in which spectra of thousands of M31 halo 
giants have been obtained (e.g. Chapman et al. 2006, Ibata 
et al. 2007) . Hypervelocity red giant stars could be detected 
in similar future surveys of stars near M31. 

There may be other sites aside from the MW and M31 
where Local Group HVSs are produced. The star discovered 
by Edelmann et al. (2005) was recently studied by Bonanos 
& Lopez-Morales (2007); the metallicity of this star is consis- 
tent with its suspected origin in the Large Magellanic Cloud. 
Due to this result, it seems likely that production of HVSs 
also takes place in the larger Local Group dwarf galaxies. 
In fact, an HVS production rate calculation for M32 by Lu 
et al. (2007) gave a rate comparable to the MW rate and 
higher than the M31 rate (though this calculation was only 
for the IBH mechanism). It is thus possible that in addition 
to the many low-mass M31 HVSs we predict to be located 
within the MW halo, there is also a significant population of 
HVSs from dwarf galaxies. This would increase the amount 
of information available from HVSs, though ambiguities in 
the origin of HVSs may make the analysis of the extragalac- 
tic HVS population more difficult. In some cases the origin 
of an HVSs could be determined from the speed and di- 
rection of its motion, but it may be difficult to distinguish, 
for instance, M31 and M32 HVSs. Of course, in order to 
make reliable predictions about HVSs from dwarf galaxies 
and determine whether they reach the MW halo in signif- 
icant numbers, detailed calculations of dwarf galaxy HVS 
production rates and distributions must be performed. We 
hope to investigate dwarf galaxy HVSs in future work. 

The use of extragalactic HVSs as probes of the Local 
Group mass distribution is certainly one of their most in- 
teresting applications. An investigation of what constraints 
these HVSs could impose on the mass distribution would be 
beyond the scope of this paper, but we hope to study this 
problem in the future. 



7 CONCLUSIONS 

In this paper we have investigated the distributions and ve- 
locities of the hypervelocity stars from the M31 Galaxy by 
integrating a large number of star trajectories within the Lo- 
cal Group. Below we summarise the main conclusions from 
our calculations. It must be noted that, due to the uncertain- 
ties in the numbers of HVSs ejected and our simple Local 
Group model, the numerical values we give are only order- 
of-magnitude estimates. 

• The differences between the galactic centres of M31 and 
the MW imply that the BHC model for HVS production 
would not be effective in M31. The TB mechanism and the 
IBH mechanism could produce large numbers of HVSs in 
M31. 

• Most HVSs ejected from M31 escape the local group, 
though some stars with lower velocities are gravitationally 
drawn towards the MW. The HVS trajectories show chaotic 
behaviour in some regions of phase space. 

• While high mass stars can only be found close to M31, 
there are large numbers of low mass (~ IMq) M31 HVSs 
on the intergalactic axis and in an arc near the MW. These 
significant overdensities close to the MW are caused by grav- 
itational focusing. We expect there to be roughly a thousand 
M31 HVSs within the virialized halo of the MW galaxy 
for the TB mechanism and a few hundred for the IBH 
mechanism. These HVSs from M31 could be distinguished 
from other stars by their large radial approach velocities 

—500 km s~^) and could thus be detected by future large 
scale astrometric surveys. 

• Our calculations of the number density of luminous hy- 
pervelocity RGB stars indicate that there are ^ 30 — 70 such 
stars in the M31 halo, of which ^ 5 should be identifiable 
and observable. 
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